Search for: All records

Abstract Radio-frequency interference (RFI) presents a significant obstacle to current radio interferometry experiments aimed at the Epoch of Reionization. RFI contamination is often several orders of magnitude brighter than the astrophysical signals of interest, necessitating highly precise identification and flagging. Although existing RFI flagging tools have achieved some success, the pervasive nature of this contamination leads to the rejection of excessive data volumes. In this work, we present a way to estimate an RFI emitter’s altitude using near-field corrections. Being able to obtain the precise location of such an emitter could shift the strategy from merely flagging to subtracting or peeling the RFI, allowing us to preserve a higher fraction of usable data. We conduct a preliminary study using a two-minute observation from the Murchison-Widefield Array (MWA) in which an unknown object briefly crosses the field of view, reflecting RFI signals into the array. By applying near-field corrections that bring the object into focus, we are able to estimate its approximate altitude and speed to be$$11.7$$km and 792 km/h, respectively. This allows us to confidently conclude that the object in question is in fact an airplane. We further validate our technique through the analysis of two additional RFI-containing MWA observations, where we are consistently able to identify airplanes as the source of the interference. 

ABSTRACT We present an analysis of Epoch of Reionization (EoR) data from Phase II of the Murchison Widefield Array using the simpleds delay spectrum pipeline. Prior work analysed the same observations using the fhd/εppsilon imaging pipeline, and so the present analysis represents the first time that both principal types of 21 cm cosmology power spectrum estimation approaches have been applied to the same data set. Our limits on the 21 cm power spectrum amplitude span a range in k space of $$|k| \lt 1 \, h_{100}\, {\rm Mpc}^{-1}$$ with a lowest measurement of Δ2(k) ≤ 4.58 × 103 mK2 at $$k = 0.190\, h_{100}\, \rm {Mpc}^{-1}$$ and z = 7.14. In order to achieve these limits, we need to mitigate a previously unidentified common mode systematic in the data set. If not accounted for, this systematic introduces an overall negative bias that can make foreground contaminated measurements appear as stringent, noise-limited constraints on the 21 cm signal amplitude. The identification of this systematic highlights the risk in modelling systematics as positive-definite contributions to the power spectrum and in ‘conservatively’ interpreting all measurements as upper limits. 

ABSTRACT We present a comprehensive simulation-based study of the bayeseor code for 21 cm power spectrum recovery when analytically marginalizing over foreground parameters. To account for covariance between the 21 cm signal and contaminating foreground emission, bayeseor jointly constructs models for both signals within a Bayesian framework. Due to computational constraints, the forward model is constructed using a restricted field of view (FoV) in the image domain. When the only Epoch of Reionization contaminants are noise and foregrounds, we demonstrate that bayeseor can accurately recover the 21 cm power spectrum when the component of sky emission outside this forward-modelled region is downweighted by the beam at the level of the dynamic range between the foreground and 21 cm signals. However, when all-sky foreground emission is included along with a realistic instrument primary beam with sidelobes above this threshold extending to the horizon, the recovered power spectrum is contaminated by unmodelled sky emission outside the restricted FoV model. Expanding the combined cosmological and foreground model to cover the whole sky is computationally prohibitive. To address this, we present a modified version of bayeseor that allows for an all-sky foreground model, while the modelled 21 cm signal remains only within the primary FoV of the telescope. With this modification, it will be feasible to run an all-sky bayeseor analysis on a sizeable compute cluster. We also discuss several future directions for further reducing the need to model all-sky foregrounds, including wide-field foreground subtraction, an image-domain likelihood utilizing a tapering function, and instrument primary beam design. 

Abstract We present deep upper limits from the 2014 Murchison Widefield Array Phase I observing season, with a particular emphasis on identifying the spectral fingerprints of extremely faint radio frequency interference (RFI) contamination in the 21 cm power spectra (PS). After meticulous RFI excision involving a combination of theSSINSRFI flagger and a series of PS-based jackknife tests, our lowest upper limit on the Epoch of Reionization (EoR) 21 cm PS signal is Δ²≤ 1.61 × 10⁴mK²atk= 0.258h Mpc⁻¹at a redshift of 7.1 using 14.7 hr of data. By leveraging our understanding of how even fainter RFI is likely to contaminate the EoR PS, we are able to identify ultrafaint RFI signals in the cylindrical PS. Surprisingly this signature is most obvious in PS formed with less than 1 hr of data, but is potentially subdominant to other systematics in multiple-hour integrations. Since the total RFI budget in a PS detection is quite strict, this nontrivial integration behavior suggests a need to more realistically model coherently integrated ultrafaint RFI in PS measurements so that its potential contribution to a future detection can be diagnosed. 

ABSTRACT In a companion paper, we presented bayescal, a mathematical formalism for mitigating sky-model incompleteness in interferometric calibration. In this paper, we demonstrate the use of bayescal to calibrate the degenerate gain parameters of full-Stokes simulated observations with a HERA-like hexagonal close-packed redundant array, for three assumed levels of completeness of the a priori known component of the calibration sky model. We compare the bayescal calibration solutions to those recovered by calibrating the degenerate gain parameters with only the a priori known component of the calibration sky model both with and without imposing physically motivated priors on the gain amplitude solutions and for two choices of baseline length range over which to calibrate. We find that bayescal provides calibration solutions with up to 4 orders of magnitude lower power in spurious gain amplitude fluctuations than the calibration solutions derived for the same data set with the alternate approaches, and between ∼107 and ∼1010 times smaller than in the mean degenerate gain amplitude, on the full range of spectral scales accessible in the data. Additionally, we find that in the scenarios modelled only bayescal has sufficiently high fidelity calibration solutions for unbiased recovery of the 21-cm power spectrum on large spectral scales (k∥ ≲ 0.15 hMpc−1). In all other cases, in the completeness regimes studied, those scales are contaminated. 

ABSTRACT High-fidelity radio interferometric data calibration that minimizes spurious spectral structure in the calibrated data is essential in astrophysical applications, such as 21 cm cosmology, which rely on knowledge of the relative spectral smoothness of distinct astrophysical emission components to extract the signal of interest. Existing approaches to radio interferometric calibration have been shown to impart spurious spectral structure to the calibrated data if the sky model used to calibrate the data is incomplete. In this paper, we introduce BayesCal: a novel solution to the sky-model incompleteness problem in interferometric calibration, designed to enable high-fidelity data calibration. The BayesCal data model supplements the a priori known component of the forward model of the sky with a statistical model for the missing and uncertain flux contribution to the data, constrained by a prior on the power in the model. We demonstrate how the parameters of this model can be marginalized out analytically, reducing the dimensionality of the parameter space to be sampled from and allowing one to sample directly from the posterior probability distribution of the calibration parameters. Additionally, we show how physically motivated priors derived from theoretical and measurement-based constraints on the spectral smoothness of the instrumental gains can be used to constrain the calibration solutions. In a companion paper, we apply this algorithm to simulated observations with a HERA-like array and demonstrate that it enables up to four orders of magnitude suppression of power in spurious spectral fluctuations relative to standard calibration approaches. 

ABSTRACT The 21 cm hyperfine transition of neutral hydrogen offers a promising probe of the large-scale structure of the universe before and during the Epoch of Reionization (EoR), when the first ionizing sources formed. Bright radio emission from foreground sources remains the biggest obstacle to detecting the faint 21 cm signal. However, the expected smoothness of foreground power leaves a clean window in Fourier space where the EoR signal can potentially be seen over thermal noise. Though the boundary of this window is well defined in principle, spectral structure in foreground sources, instrumental chromaticity, and choice of spectral weighting in analysis all affect how much foreground power spills over into the EoR window. In this paper, we run a suite of numerical simulations of wide-field visibility measurements, with a variety of diffuse foreground models and instrument configurations, and measure the extent of contaminated Fourier modes in the EoR window using a delay-transform approach to estimate power spectra. We also test these effects with a model of the Hydrogen Epoch of Reionization Array (HERA) antenna beam generated from electromagnetic simulations, to take into account further chromatic effects in the real instrument. We find that foreground power spillover is dominated by the so-called pitchfork effect, in which diffuse foreground power is brightened near the horizon due to the shortening of baselines. As a result, the extent of contaminated modes in the EoR window is largely constant over time, except when the Galaxy is near the pointing centre. 

Abstract Precise instrumental calibration is of crucial importance to 21-cm cosmology experiments. The Murchison Widefield Array’s (MWA) Phase II compact configuration offers us opportunities for both redundant calibration and sky-based calibration algorithms; using the two in tandem is a potential approach to mitigate calibration errors caused by inaccurate sky models. The MWA Epoch of Reionization (EoR) experiment targets three patches of the sky (dubbed EoR0, EoR1, and EoR2) with deep observations. Previous work in Li et al. (2018) and (2019) studied the effect of tandem calibration on the EoR0 field and found that it yielded no significant improvement in the power spectrum (PS) over sky-based calibration alone. In this work, we apply similar techniques to the EoR1 field and find a distinct result: the improvements in the PS from tandem calibration are significant. To understand this result, we analyse both the calibration solutions themselves and the effects on the PS over three nights of EoR1 observations. We conclude that the presence of the bright radio galaxy Fornax A in EoR1 degrades the performance of sky-based calibration, which in turn enables redundant calibration to have a larger impact. These results suggest that redundant calibration can indeed mitigate some level of model incompleteness error.